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Abstract 

We compute the shear viscosity and the electrical conductivity in gauge 
theories with massive fermions at leading order in the large Nf expansion. 
The calculation is organized using the 1 /Nf expansion of the 2PI effective 
action to next-to-leading order. We show explicitly that the calculation is gauge 
fixing independent and consistent with the Ward identity. We find that these 
transport coefficients depend in a nontrivial manner on the coupling constant 
and fermion mass. For large mass, both the shear viscosity and the electrical 
conductivity go to zero. 
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1 Introduction 


Transport coefficients in relativistic gauge theories have been discussed in a number 
of papers in the past few years m las s mm mi- The motivation comes from 
possible applications in heavy ion physics and the early universe, as well as from the¬ 
oretical interest. However, the attention has mostly been focused on ultrarelativistic 
theories, where the scale is set exclusively by the temperature. In this paper we un¬ 
dertake the computation of transport coefficients in gauge theories at temperatures 
where the fermion mass cannot be neglected. We carry out this study in the large 
Nf limit of QED and QCD, where a complete leading order calculation is possible. 
For massless fermions transport coefficients have been computed in large Nf gauge 
theories in Ref. 0, using kinetic theory. A study of thermodynamic properties of 
gauge theories in the large Nf limit can be found in Ref. PI 

We perform a diagrammatic calculation, organized using the 1 /Nf expansion of the 
two-particle irreducible (2PI) effective action to next-to-leading order (NLO). The 2PI 
effective action is a useful tool in studying the nonequilibrium dynamics of quantum 
holds m ■ In actual applications, the 2PI effective action is truncated at some order in 
a chosen expansion parameter. In Ref. []2j it was shown for a number of theories that 
the lowest nontrivial truncations correctly determine transport coefficients at leading 
(logarithmic) order in the expansion parameter. Here we show explicitly that the 
lowest order nontrivial truncation of the 2PI effective action in the 1 /Nf expansion 
provides all the required ingredients to successfully compute the shear viscosity and 
the electrical conductivity. When considering gauge theories and effective actions, 
care is required with respect to gauge invariance and Ward identities [T2J • We show 
that despite the nontrivial resummation of diagrams carried out, the method provides 
a gauge fixing independent result and is consistent with the Ward identity. This 
provides an explicit example of a nontrivial quantity for which potential non gauge 
invariant contributions in a fully self-consistent calculation would be suppressed by 
powers of the expansion parameter. 

The paper is organized as follows. In Section [21 we formulate the 2PI effective 
action to NLO in large Nf QED. We obtain the integral equation relevant for the 
calculation of the shear viscosity and the electrical conductivity and discuss power¬ 
counting in the 1 /Nf expansion. We show that a typical diagram that contributes 
to the shear viscosity at leading order in the large Nf expansion is as shown in 
Fig. |H Plasma effects relevant for transport coefficients are studied in Sectional In 
the following Section, we explicitly work out the integral equation relevant for the 
shear viscosity and write it in a form convenient for a variational treatment. In Sec¬ 
tion El this analysis is repeated for the electrical conductivity and the Ward identity 
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Figure 1: Typical skeleton ladder diagram that contributes to the shear viscosity in 
large Nf QCD at leading order in the 1 /Nf expansion. 


is explicitly checked. We generalize the discussion from QED to large Nf QCD in 
Section |H1 The numerical analysis and our results are presented in Section [TJ The 
final Section is devoted to the conclusions. In Appendix El we derive a set of integral 
equations from the 2PI effective action which are employed in the main text. Finally, 
Appendix El contains parametric estimates in the leading logarithmic approximation, 
for both massless and very heavy fermions. 

A short summary of these results has appeared in Ref. m Part of the analysis 
is very similar to the study of the shear viscosity in the O(N) model in the large N 
limit ns- When possible, we will refer to that paper for further details. 


2 2PI-1/AT expansion 


Since the structure of QED and QCD is similar in the large Nf limit, we use QED 
in the following for the purpose of discussion. Color factors will be introduced later. 
The action for Nf identical fermion fields i/ a (a = 1,... ,Nf) then reads 1 


S 



l -F^ + i, a {ip 


m) i\) c 


+ Agf + S, 


gh; 


( 1 ) 


with 


P = 


D U ~ dn + 


ie 




An, 


( 2 ) 


and we use the notation 



(3) 


where C is a contour in the complex-time plane. Note that we have rescaled the 
coupling constant with ~Nf , so that in the large Nf limit Nf goes to infinity while e 

1 We use = diag(+, —, —, —), so that P 2 = — p 2 , p = |p|. The y-matrices obey {7C7Q = 

2 g^ u . Traces over Dirac indices are indicated with tr. 
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Figure 2: NLO contribution to the 2PI effective action in the 1 /Nj expansion. 

remains finite (after renormalization). To fix the gauge we use a general linear gauge 
fixing condition, 

Sg = - l (4) 

Below we specialize to the generalized Coulomb gauge: /o = 0, /, = <9*. The ghost 
part is not needed explicitly. 

The 2PI effective action is an effective action for the contour-ordered two-point 
functions 

D^(x,y) = {T c Ap(x)A„(y)), S ab (x, y) = (T c ijj a (x)jp b (y)), (5) 

and can be written as |T3 

F[S,D] = ^Trln D~ l + ^Tr D^{D - D 0 ) 

—zTr In S '” 1 — iTr Sq 1 (S — S 0 ) + r 2 [S', D] + ghosts, (6) 

where Dq 1 and Sq 1 are the free inverse propagators. The 2PI effective action frame¬ 
work automatically entails the existence of a set of coupled integral equations for the 
various 4-point functions. These integral equations contain the relevant physics for 
the calculation of some transport coefficients in a number of theories JO]- In Ap¬ 
pendix El we briefly describe how to obtain the relevant set in the theory we study 
here. 

The lowest order contribution to r 2 [>S', D\ appears at next-to-leading order (NLO) 
in the large Nf expansion (see Fig. II 

r ^ LO [S,D] =[ trYSab{x,y)YSba{y,x)D flu (x,y). (7) 

Zly f Jxy 

We specialize to the completely symmetric case and write S ab = 5 ab S, T, ab = 5 ab £• 
The resulting self energies (see Fig. ED are then 

W"{x,y)=eHx 1 *S(x,y)'fS{y,x), 

Z(x,y) = —^-'fS(x,y)'fD lu/ (x,y). 
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( 8 ) 

(9) 




Figure 3: Fermion and gauge boson self energy. 


They depend on the full propagators, determined by the Dyson equations 

zr 1 = d - 1 - n, s ' -1 = Sq 1 - s. (10) 

The set of integral equations for the 4-point functions (see Eq. (11441) in Ap¬ 
pendix 0, up to this order in the large Nf expansion, is shown in Fig. 01 These 
coupled equations sum all the diagrams that are required to obtain the shear viscos¬ 
ity and electrical conductivity at leading order in the 1 /Nf expansion |12j . This can 
be argued as follows. Kubo formulas relate these transport coefficients to the slope 
of current-current spectral functions at vanishing frequency 


Id , n . 


Id , o . 

° . 0 ) 


q°=0 6 dq° JJ ' q°= o' 

where the spectral functions are defined as 

Pmr{x -y) = <[t T ij (x),7l ij (y)}), Pjj(x ~ y) = ([/(x),/(?/)])• 

Here 7 is the traceless part of the spatial energy-momentum tensor, 


Kij(x) ^ijFk»F k » i'lpc 


(ii) 


( 12 ) 


( jiEi _ 2l£i _ 1 5 in „D^ V«, (13) 


and j l (x) = qf'tjj a (x)Y'ijj a (x) is the electromagnetic current, with qf the charge of the 
fermion. 2 The correlators in Kubo formulas are computed in thermal equilibrium, so 
from now on we specialize to the Matsubara contour and work in momentum space. 

The correlators in the Kubo formulas are required in a specific kinematic config¬ 
uration which, as is well known, causes them to suffer from so-called pinching poles. 
These pinching poles are screened by the imaginary part of self energy, leading to the 
appearance of a factor inversely proportional to this imaginary part. This modifies 

2 This is the charge with which the fermions couple to the external operator; we prefer to distin¬ 
guish it from the coupling between the gauge bosons and the fermions in the ladder diagrams. In 
QED it is also rescaled, so that qf = e/y/Nf, while in QCD it is not. 
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Figure 4: Integral equations for the 4-point, functions at NLO in the 1 /Nf expansion. 


the naive power counting scheme. The fermionic one loop diagram, which contributes 
to both transport coefficients, is naively of order Nf, due to the Nf identical fermion 
fields that run in the loop. Because of the pinching poles, this is enhanced by the 
inverse thermal width (given by the imaginary part of self energy) which is of order 
1 /Nf, as we show below. We find therefore that the conductivity and the shear vis¬ 
cosity are proportional to Nj in the large Nf limit (apart from the external charges in 
the case of the electrical conductivity). Adding a vertical photon line to the one-loop 
fermion diagram gives a contribution that is of the same order; the factor of 1 /Nf 
from the added vertices is compensated by a new pair of propagators with pinching 
poles and a new inverse factor of the thermal width. This remains true when adding 
any number of vertical photon lines; therefore all these diagrams have to be summed. 
A contribution of the same order is also obtained when considering a box rung with 
horizontal photon lines and vertical fermion lines (see Fig. n. In this case, a new 
pair of propagators with pinching poles along with a new closed fermion loop com¬ 
pensates for the additional four coupling vertices. Note that there are two ways a 
box rung can be added, depending on the orientation of the fermion lines. Again, a 
diagram with any number of box rungs contributes at leading order as well. These 
kind of diagrams are precisely those which are summed by the integral equations for 
the fermionic 4-point function we obtained from the 2PI effective action. 

In the case of the shear viscosity, the external operator also couples to two gauge 
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Figure 5: Rungs in the integral equation for the fermion 4-point function. 


boson fields. It is therefore necessary to consider diagrams with gauge bosons on the 
side rails. Again, the corresponding imaginary part of the gauge boson self energy 
has to be included in the side rails propagators to avoid pinching poles. If the gauge 
boson is an onshell stable excitation, its self energy in Fig. [3 yields a thermal width 
only when at least one of the fermion lines in the diagram is dressed, a contribution 
of order 1/Nf (see Eq. As a result the gauge boson thermal width is of order 

1/Nf, similar as the fermion thermal width. On the other hand, if the gauge boson 
is an unstable excitation, no fermion lines need to be dressed in the gauge boson 
self energy to get a non-vanishing imaginary part, which is therefore of order Nf ; 
in this case pinching poles do not lead to a further enhancement. However, in both 
cases there is only one gauge boson compared to Nf fermion fields. Therefore these 
diagrams are subleading in the 1/Nf expansion. For the shear viscosity we finally 
also have to consider diagrams where one external operator couples to two gauge 
boson fields and the other one to two fermion fields, and there is at least one fermion 
rung. In this case there are pinching poles from the pair of gauge boson propagators 
and from the pair of fermion propagators. If the gauge boson is an onshell stable 
excitation, we find two powers of Nf from the pinching poles, one power of Nj from 
the closed fermion loop and a power of 1/Nf from the coupling vertices. However, 
due to kinematics this diagram gives a nonzero contribution only when the spectral 
density of the fermionic rung is offshell, which introduces a further power of 1/Nf 
and makes the contribution from this diagram subleading. If the gauge boson is an 
unstable excitation, the fermionic rung can be onshell but we find only one power of 
Nf from the pinching poles and the diagram is subleading as well. We conclude that 
diagrams where one or both of the external operators couple to gauge bosons can be 
neglected at leading order in the large Nf expansion. 

It is therefore not necessary to consider the full set of integral equations in Fig. 01 
only the closed integral equation for the 4-point, function where all external legs are 
fermionic is required. In this respect, the large Nf computation is slightly easier 
than the leading-log calculation in the weak-coupling limit, where two coupled inte¬ 
gral equations for the fermion and the gauge boson contributions have to be solved 
simultaneously i3j. Instead it is very similar to the analysis in the O(N) model in the 
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large N limit, with the gauge boson and the bubble chain playing a similar role DU- 
The individual kernels in the integral equations in Fig. 0 are obtained by cutting 
one line in the self energies and read 

A a b ;cd(R, P) = -^-S ad 5 bc YD^(R - P)r, 
iv / 

P) = - P)Y + YS(R - P)-f] . (14) 

iV / 

A „„AR. P ) = h“S(p - R)Y + YS(P - R)-f], 

iv / 

where R is the momentum that enters and leaves on the left and P enters and leaves 
on the right. To obtain a closed integral equation for the fermionic 4-point function, 
the third equation in Fig. 0 is substituted into the first one, leading to 


r K ) = W*. K) - f p aP) 5(P)5(P)r^(p, k), (is) 

with the effective kernel 


A ab;cd{R, P) — A a b;cd{R, P ) + “ ^ A L) D ua (L) D pl j,(L) A pcr;c d(L, P ) 

l\f 


(16) 


+^ ab <5 cd [fW(P - L)^ + rS(R + L)Y] D va {L)D p ^L)YS{P - L) 7 ° 


We use the notation 


t r T E/. 

7p „ d P 


d 3 p 

(2 


(17) 


where the sum runs over the corresponding Matsubara frequencies. To carry out the 
frequency sums, it is convenient to introduce a 3-point effective vertex r ab (P + Q, P) 
as 

r„ t (p + q,p) = C( P ) ~rf, R S(R + Q)rh(r)s(R)r%jR. p- q), (is) 

where F° fe (p) is the bare coupling between the fermion fields and the external operator 
under consideration, and Q is the momentum entering the operator insertion. This 





Figure 6: Integral equation for the full 3-point function. 


yields the final integral equation we work on in the remainder of this paper (see Fig. m 

r ab (P + Q,P) = r° 6 (p) - ^ S(R + Q)T cd (R + Q, R)S(R)A cd . ab (R, P; Q), (19) 

with the kernel 

2 4 

A ab; Cd (R, P; Q) = ~8 ad 8 bc YD^{R - P)Y + ^8 ab 8 cd (20) 

x ^ [YS(R - L)r + 7 »S(R + L + Q)Y] D^L + Q)D m/ (L)YS{P - L ) 7 CT , 

A typical skeleton diagram that contributes to the shear viscosity in large Nf gauge 
theories is depicted in Fig. [Q Throughout the paper we drop subleading powers of Nf. 


3 Quasiparticles 

In this section we study the effects of the medium on the propagation of both fermions 
and gauge bosons at this order in the 1 /Nf expansion. In particular, we discuss the 
fermionic thermal width and compute the full gauge boson self energy required at 
this order in the Nf limit. 


3.1 Gauge boson 


We choose to work with the gauge boson propagator in the generalized Coulomb 
gauge, so that it reads 


pi< pi' 

D^{P) = PTAt(P) + ^V°Al(P) + 

with transverse and longitudinal propagators 

A (p) = _1_ r° duj p T {uj, p) 

T ujl + p 2 + n T (P) J 2vr iu n — uj ’ 

A (p) = - 1 = _ r duj p L (uj, p) 

L p 2 + II L (P) p 2 /-«> 2vr iu n - cj ' 


( 21 ) 


( 22 ) 

(23) 
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In this gauge the gauge boson spectral function is independent of £, 


The self energy 

is decomposed as 


lT(P) = Pf-p T (P) + s'Wd’)- 

W(P) =e 2 X tr-r “S(P + K)YS(K), 


K 


P 2 


n ^( p ) = p^u t (p) + —p^Ul(p), 


p* 


with the usual projectors 


Pft = S ij - ftp?, P£° = Pft = 0, 
The transverse and longitudinal self energies 

n L = -n 00 , n T = 


Pft = g^ 

are then 


pn pv 

P 2 


+ Ptft. 


nc + 


p 2 

ft 2 


-n 1 


00 


(24) 

(25) 

(26) 

(27) 

(28) 


Since we drop subleading corrections in the 1 jNf expansion, and pinching poles are 
not an issue here, the fermionic propagators in Eq. m can be taken at leading order, 
i.e. free ones. 

We need to compute both II 00 and lift We split the self energy into vacuum and 
thermal parts 

IT- = IC, + rc (29) 

where the vacuum contribution has the usual form 


n £(P) = ( P 2 g““ - w) n TO (p), 


(30) 


with 

IWP) = () + ln4 >r - 7 b) + n L ( -P)- (31) 

We used dimensional regularization in 3 — 2e dimensions and eo is the bare coupling 
constant. In order to carry out the renormalization, 3 we introduce the dimensionless 
running coupling constant in the MS scheme 


1 _ h 2e , 1 

e 2 (/i) e 2 12vr 2 


+ ln47r — 7 e 


3 For renormalization of 2PI effective actions beyond what is needed here, see Ref. ca 


(32) 
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The running coupling constant obeys the usual renormalization group equation with 
/3(e 2 ) = e 4 /(67r 2 ). Renormalization is now straightforward, 

/i~ 2e e 2 A^ = e 2 (/i)A T/L , (33) 


with the renormalized propagators, 

1 


At 1 — 


-p 2 (i + n 


/ i 

vac ) 


+ n^ h ’ 


Ar, — 


— 1 


P 2 (l 


n / i 

vac J 


+ m h 


(34) 


We note here that the product e 2 (p)/S .t / l is renormalization group invariant. As is 
well known, the theory has a Landau pole at the scale A l = // e 67r2 / e2 (/''). The largest 
scale in the problem, either the temperature or the mass, has to be reasonably well 
below the Landau scale. This imposes a restriction on the allowed values of the 
coupling constant. Although the results are renormalization group invariant, in order 
to present them numerically we have to choose a scale. To facilitate a comparison 
between our results and the ones obtained for massless fermions in kinetic theory [2j, 
we take p = p dr = 7 re -7B T, the dimensional reduction value for massless fermions. 

The real part of the finite piece at zero temperature reads 


RenL(P) = 


12ir 2 


| [0(P 2 - 4m 2 ) + 0(-P 2 )] (l + P(P) 


In 


1 - (3{P) 


-20(4 m 2 - P 2 )0(P 2 ) 1 + 


2 m 2 
~P* 


B (P) arctan 


1 + P(P) 

1 


P(P) 


4m 2 

”W~ 


, m 2 5 

V + 5 


where 


P(P) = \ 1- 


4m 2 


P(P) = 


4m 2 


- 1. 


p 2 ’ “C / U p2 

We now consider the thermal piece. For the real parts we find 


Ren !: RU p> = 

Os 

to 

7T 2 

Ren “, h (p) = 

2e 2 

7T 2 


f°° k 2 

/ dk—npiuj k) 
'o ^k 

roc ^ 2 

/ dk — ufPP) 

' 0 ^k 


P 2 + 2 m 2 , 

1 H-—-In 


1 + 


8 kp 

P 2 + 4o^ 2 

8 kp 


(k + p+)(k + p ) 


In 


(k - p + )(k - p_) 
(k +p + )(k +p 


P°u k 
2 pk 


In 


(k~P+)(k~P-) 
(P 2 — 2p°u k — 2 pk)(P 2 + 2p n w k + 2 pk) 


(P 2 — 2pV k + 2 pk){P 2 + 2pV k — 2 pk) 


(35) 


(36) 


■ (37) 


(38) 


11 

































Figure 7: Transverse and longitudinal spectral functions pt,l(^,p) for p/T = 2, 
m/T = 2 and g e s = e = 2. For these parameters, the transverse gauge boson is a 
stable quasiparticle at Ut/T = 2.05, indicated with the vertical line. 


where p± = ~[p±p°/3(P)]. The remaining one-dimensional integrals can be done nu¬ 
merically. The imaginary parts, including the vacuum contribution, can be computed 
explicitly and read 4 


Im ITT (P) 

PR, th v ' 


: — (P 2 + 2m 2 )(©(P 2 
47T [ 

+0(-P 2 ) 


4m 2 ^ 

2 T 
P 


2 T 1 4- p-P+/ T I 

m + - In + % 

p 1 + e p -/ 1 

1 + e ~P+/ T ] 

In 4 


1 + eP-/ T 


( QQl 


Imn“ th (P) 


e 2 T 2 


/ 7 ~ 


+/3(P) [Li 2 (-e^+/ T ) + Li 2 (-e- p_ -/ T )] + y [Li 3 (-e-^/ T ) - Li 3 (-e^-/ T )] | 

4 The arguments of the polylogarithmic functions Li n (z) are chosen such that they lie between 
— 1 and 0 for positive p°. 
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Figure 8: As in Fig. 0 with m/T = 0.1. In this case the transverse gauge boson is a 
resonance in the continuum at ujt/T = 2.13. 


e 2 T 2 


7r 


m 2 p 




1 + e~ p +/ T 


P 2 T 
2 T 
V 


+ /3{P) [U 2 {-e- p +' T ) - U 2 {-e p ~/ T )] 


1 + e p -l T 
[Li 3 (—e _p+/T ) - Li 3 (-e p - /T )] 


(40) 


where p± = \ \p° ± p/3(P)\. 

The resulting transverse and longitudinal spectral densities are shown in Figs. [7| 
and|H The propagating transverse and longitudinal modes, ujt/l(p), a re determined 
from the poles of the corresponding propagators in Eq. (ED- For the specific param¬ 
eters chosen in Figs. 0and|Hl the transverse gauge boson is stable in the first case 
and a resonance in the second. In both cases the longitudinal gauge boson does not 
propagate. 
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3.2 Fermion 

The fermionic propagator is 

S(P) = T 


TUnT 


7 • p — m — E(P) 


du p F (uj,p) 
2n iCur, — u ' 


(41) 


where u n = (2 n + 1 ) 7 rT (n E Z) is a fermionic Matsubara frequency and Pf(w, p) the 
spectral density of the fermion. The self energy can be decomposed as 


E(P) = 7 °S°(P) + 7 • p S S (P) + £ m (P), 
such that the retarded electron propagator reads 

7° |P° - Z° R (P)} - 7 • P b> + Sit/ 5 )] + [m + £g(P)] 


S«(P) = 


[po _ E» r (P )] 2 - (p + 57-P )] 2 - [m + £?(■?)]= 


(42) 


(43) 


The poles of the retarded propagator at p° = E p — iT p /2 determine the properties of 
the quasiparticle excitations of the system. The self energy is 


s(p) = - Wf y R rs(R)rD P „(R-p), 


(44) 


and we find 


E p 


e> — 

N f 


T P = 


ImS^(P) 


p u 


p°==bajp 


+ ol 7 


(45) 


with cn p = \Jp 2 + m 2 and 


=S(P) = 2 b°S° K (P) + pS p R (P) +mE“(P)] = -trE R (P)(f + m). 

The retarded and advanced propagators then simplify to 

f + m 


^/a(P) = 


P 2 — m 2 — Tin E- / 4 (P) ’ 
and the corresponding spectral density is 

—2ImE^(P) 


(46) 


(47) 


Pf(P) = (f+ m) 


+ [Im£g(f>)]' 


(48) 
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In the large Nf limit, this reduces to the spectral density of a free fermion, 


Pf(P ) = + m) p(P) = iff + m) 2 tt sgn(p°)S(P 2 — m 2 ). (49) 


Nonetheless, whenever a pair of propagators with pinching poles is present, the large 
Nf limit of the product of a retarded and an advanced propagator is 

S r (P)YSa(P) =[p + m) Y (f + m) + C>(jV»). (50) 

It remains to give the explicit expression for the thermal width. Combining 
Eqs. da BSD we find, after doing the Matsubara frequency sum using spectral repre¬ 
sentations, performing the analytic continuation and taking the trace, 


2e 


Imtr£ R (P)(f + m) = — — / [n F (r°) + n B (r° - p 0 )] p(R) 


NfjR 


x \ 2pr(R-P) 


r°p° — m A — r • k p • k 


+ Pl{R~ P ) [r°P° + m 2 + r ■ p] | , (51) 


where k = r — p. Since only the spectral density of the gauge boson propagator 
contributes to the thermal width, this is explicitly independent of the gauge fixing 
parameter. 

We proceed by introducing k — |r — p| as 


1 



r “Pl) 



k 

dk —5(cos 6 pr 
rp 


Zpr ) 5 


(52) 


where cos 0 pr = p • r is the cosine of the angle between p and r and 


Zpr 


r 2 + p 2 — k 2 


2 rp 


(53) 


The integral over r° is performed using p(R) and the one over 6 pr using the delta 
function introduced above. After that, the final result for the thermal width reads 


r p 


2 poo pr+p 

- 5 - / dr— dkk 

8N f ir 2 pu p J o ay J\ r _ p \ 

x | [n F (u r ) + n B (u r + (Up)] [2 c^p T (uj r + cu p , k) - c^p L (u T + cn p , k)] 

- [n F (uJr) + n B (uJr - o; p )] [2c^pr(^r - Wp, k) - clp L (uj r - (Up, k)] i, (54) 
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Figure 9: Effective one-loop diagram contributing to the shear viscosity. 


with 


, pr 2 

c T = ±bj r uj p + ^ (r - pz pr )(p - rz pr ) + m 


c± = =Fu; r u;p + prz pr + m 2 . (55) 


It should be noted that the thermal width as it is written here is not defined, due 
to the divergent contribution from soft quasistatic transverse gauge bosons HBJEH!. 
However, in the application to transport coefficients this contribution cancels against 
part of the ladder diagrams. This has been analyzed in detail in the weak coupling 
limit in Refs. m 


4 Shear viscosity 

We are now ready to compute the shear viscosity. It is obtained from the effective 
one-loop diagram (see Fig. H 

GWQ) = -5jtTS(P + Q)T i *(P + Q,P)S(P)T%(p), (56) 

where Q = (iui q , 0 ) is the momentum that enters from the left. The coupling between 
the external operator and two fermionic fields is Tg (p) = p T 13 py, p], with 

T lJ [7, P] = \ {YP J + 7 J P l ) ~ ' P- (57) 

In order to carry out the sum over Matsubara frequencies, we use the method de¬ 
scribed in Ref. |4]. From the analytic properties of the propagators and the structure 
of the integral equation, it follows that the effective vertex Y l ' J (P + Q,P ) has cuts 
along Im(p°) = 0 and Im(p° + q°) = 0. As is usually the case 4,^], To| . hr the pinching 
pole limit only one particular analytic continuation of the effective vertex is required. 
At leading order in the l/Nf expansion we arrive at 

lim p 7r7r (g°, 0) = -2 q°N f [ n' F (p°)Retr S R (P)T l:, (p 0 + i0,p° - p) S A (P)T% (p). 

9 0 -* o J P 

(58) 
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Using now Eq. m for the product of the retarded and the advanced Green functions, 
the viscosity reads 

V = + m)T lJ (p° + iO,p° - + m)T l 0 \p). (59) 

Due to the pinching poles, the momentum in the loop is forced onshell: p° = iutp. 
We may therefore decompose If* + m into spinors according to 


m 


p°=Up 


^u A (p)u A (p), 


If + m 


A=1 


p°=-u p 


-J]w(-p)w(-p), (60) 


A=1 


and associate the spinors with the effective vertex Y l ' J as 
Re h A (p)(a;p, cUp ; p)wa(p) 

= Reh A (— p)r*- ? (— o; p , —ut p ; p)u\(—p) = 2p 2 T lJ [p, p\V(p). 

An equivalent expression for the resulting scalar vertex V{p) is 


(61) 


T>(p) — —I’e i.r r’ y {p’ + i0,p° — 10; p)(f + rn}T‘°'p. p] 
o p z 


P°=±LU p 


(62) 


The normalization is such that for the bare vertex Tq (p) this yields T>(p) = 1. It is 
then straightforward to obtain 5 


4iV/ f P 2 , , \ / \ 2 N f 

= - tep 


dp —n' F (ui p )x(p), 


UJ r 


where 


xip) = 


P 2 V{p) 

r p 


(63) 


(64) 


5 A different route to arrive at these expressions is as follows. Inspection of the integral equation 
shows that in the pinching pole limit the effective vertex remains linear in the y-matrices and the 
identity and can be taken traceless. It can therefore be decomposed as 

RRP) = A(P)pT%, P] + T ij [p, p] {P(P)y • P + G(P) 7 V 0 + D(P)m} , 

with 4 independent scalar functions A, P, C and D. However, an analysis of the integral equation 
with this decomposition shows that only one linear combination of those four functions enters in the 
final scalar equation, namely 


1 


V{p) = —Re { p 2 [A(p) + B(p) + C{p)\ + to 2 [C(p) + D(p)}} 


V 


p°=±u 


as expected from the spinor decomposition. 
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4.1 Integral equation 

From the integral equation (USD, we now obtain an integral equation for the effective 
vertex T>(p). To do the frequency sums we proceed as before (see Ref. [T^] for more 
details in a similar computation). We find 

V(p) = 1+ J \ n B(r° - p°) + n F (r 0 )] ^P 2 (p ■ r)-^^-V(r)A hne (R, P) 

+ J ^ [n B (r° - p°) + np(r 0 )] ^P 2 (p • r)|^-P(r)A box i(P, P) 

+ J \n B (r° + p°) + n F (r°)\ } -P 2 {p-i)^^V(r)A hox2 (R,P), (65) 

where here and below p° = ±cu p . The three contributions arise from the line diagram 
and the two box diagrams respectively. The second Legendre polynomial P 2 (x) = 
(3x 2 — l)/2 originates from 

T lJ [r,r]T v [p,p] = jjp 2 (p-r). (66) 

The three kernels read 

Aim ,(R, P) = Yj tr h l ‘(f + m h"<f + ™)] p^(R - P), 

Atom (R,P) = 4f,j L “ r0) “ Ml ° ~ p0) l D f^ L )Dt(L) 

xtr [Yffi + m)YP f{R - L)\ tr [' Yp F {P - L)j a (]J > + m )\, 

Ab. 2 (S,P) = j^ / [n F (l a + r 0 )-nr(l°-p°)]D«(L)D?„(L) 

f J L 

xtr [ 7 "(# + m)Yp F (R + L)} tr [Ypf(P ~ L) y CT (f + m)\ . (67) 

Using Eq. (EED for the fermionic spectral functions, it is easy to see that under a 
change of variables R —» — R, the contribution from the second box diagram becomes 
identical to the first one. We can then write 

V(p) = 1 + J [n B (r° — p°) + n F (r 0 )] ^P 2 (p • r)^^-P(r)A(P, P), (68) 

with 

HR, R) = y^; tr + ™)a v (F + m )] pH r - R ) 


18 










+ W l K(/ - - n X r ° - '“)] P<- P - L ^ R - L ) D %( L )°t(L) 

f J L 

xtr [ 7 ^(Ifl + m )+ m )\ U [l p {f + m ) 7 ff (f > + m )\ ■ (69) 

In terms of the function xip) defined in Eq. (1(141) the integral equation reads 

^P r p xip) = P 2 + \ [ [nB(r° - p°) + n F (r 0 )] P 2 (p ■ r)^x(r)p(P)A(P, P), (70) 

z Jr r 

where p° = ±cn p . Upon solving it for x(p), we obtain the shear viscosity from Eq. (T77TT)l . 


4.2 Variational approach 

Since the integral equation looks prohibitively difficult to solve analytically, we pro¬ 
ceed to formulate it as a variational problem, which gives a convenient formulation 
for finding a numerical solution md 
After multiplying Eq. (EDI) with 


P ' ( ^ 

n F\ ijJ p)i 

U)p 


the integral equation can be written as 


f(p)x(p) = S(p) + 


drH^p, r)x(r) 


(71) 


(72) 


with 

4 

P(p) =p 2 n' F (uj p )r p , S(p) = —n' F (uj p ), (73) 

CJp 

and a symmetric kernel, 7i(p,r) = 77(r, p), whose explicit form is presented below. 
Since 7 i is symmetric, Eq. (d can be derived by extremizing the functional 


Q[x] = dp 

Jo 


1 9 1 f°° 

S(p)x(p) - ^(p)x 2 {p) + 2 f dr? t(p, r )x(r)x{p) 


The viscosity is then given by the extremum of the functional 

4 N f 


V = 


15n 2 


Q[X Xext]- 


(74) 


( 75 ) 


In the rest of this section we explicitly evaluate H(p,r) 
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We separately compute the single line diagram and the box diagram, TL = Ttn ne + 
7"4>ox- We start with the diagram containing the single line. Proceeding as we did in 
the calculation of the thermal width, we arrive at 6 


e (p,r) = 


e p r 

—n F (u p )- 


cr+p 


dk kP 2 (z pr ) 


8Nfir cjp to Y J\r—p\ 

x| [n F (u r ) + n B (u r + a>p)] [2c^p T (u r + u p , k ) - c£p L (cu r + u p , k)] 

+ [n F (u r ) + n B (u r - o; p )] [2c F p T (u r - u p , k) - clp L (uj r - u p , k)] 


(76) 


The coefficients crp/ L were already defined in Eq. (1551) . The contribution from the line 
diagram is £ independent for the same reason as the thermal width. Using similar 
properties of the distribution functions as discussed in Ref. (T^j , it is straightforward 
to verify that 7~fii, 1P is symmetric under exchange of p and r. 

For the contribution from the box diagrams we work out the traces and contrac¬ 
tions in Eq. (E2D- After a bit of algebra we End' 


D ?p( L ) D t( L )to -ty + rn)\ tr [ 7 p (f -If + m) 7 CT (f + m)\ 

= Af (L) A a l (L) (2tq - r°l° - r -1) (2pg - pY - p • l) 

+ 2Re[A£(L)A£(L)] (2r° - 1 °) (2p° - Z°) (p • r - p • 1 r • l) 


2A F (L)A F (L) R ■ L P ■ L - P ■ L ( r 2 - (r • l) 2 ) - R ■ L (p 2 - (p ■ 1) 


+2 (p • r — p • lr • 1 


Z (2R ■ L - L 2 ) (2P ■ L - L 2 ) 2£ (p • r - p • lr ■ l) ReA^(L) 


P 


e 2 


+2£r u p u ReA £(L) + -jjR ■ LP ■ L 


(77) 


The terms that depend on the gauge fixing parameter £ are proportional to 2 P-L — L 2 . 
These terms are accompanied by the Dirac delta function p(P — L ) (see Eq. (1551) 1. 
which for onshcll momentum P causes this factor to vanish. Since this was the last 
piece that depended on the gauge boson propagator, we hnd that the viscosity is 
gauge fixing independent, as it should obviously be. 

To proceed further, we consider the 8-dimensional integral over R and L in Eqs. 
(1701551) . The cosine of the angle between p and 1 is denoted as cos 9 p i, between r and 

6 In Ref. jT5], the first term in Eqs. (5.17) and (5.32) is written with the wrong sign. 

'In Appendix B of Ref. 0, the factor Re appears written erroneously as |A^A^|. 
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1 as cos 9 r i, and the azimuthal angle between the p, 1 plane and the r, 1 plane as qt>. 
We specialize to p° = u> p . The 8-dimensional integral can then be written as 


2t r 

(2vr) 8 





dcos 9 p i 


dcos 9 r i 


i -i 



(78) 


The integration over cos 9 p i will be performed using the delta functions in p(P — L), 
the one over cos 6 r i using p(R — L ) and that over r° using p(R). The product of the 
three spectral functions yields a set of constraints, since 


p(R)p(P-L)p(R~L) 


P°=U P 


y; <5(r° + Sia; r )(5(a;p - 1° + s 2 uj p _i)5(r 0 - 1° - s 3 a; r _i) 

Si=± 

5(a7 p + Siu; r + S2(n p _i + S3u; r _i). (79) 


Sj=± 


Out of the eight combinations, only three can contribute for kinematical reasons, 
namely those corresponding to 2 2 processes. We treat these three cases separately 

and write 


Tfbox = 'H 


(i) 

box 


+ 

+ '‘'box 


+ 77 ^ 

+ '‘'box' 


(80) 


1. (si, s 2 , s 3 ) = (—, +, —). The cosines are cos 9 p i = z p i, cos 9 r i = z rl , where 


l 2 -1 2 uj° 


Zpl 


+ 


/ 2 _ 72 

Z rl = 


u; r / 0 


2pZ ' pi ’ ~ r( 2r7 rl 

The constraints from the spectral functions can be satisfied provided 


1° > V/ 2 + 4M 2 , |/_| < p,r < |/ + |, 

where we use again the notation 

i± = 5 f ± ‘°m\, m = \p 


AM 2 


L 2 ' 


(81) 


(82) 


(83) 


The angle 0 appears both in expression (1771) and in P 2 (p ' r) in Eq. (17(1 . since 
p ■ r = sin dpi sin 9 r i cos 0 + cos 9 p i cos 9 r i. (84) 

Performing the integral over 0 yields an expression of the form 


cll |A£(L)T + c TL Re A«(L)A£(L) + c TT IA"(L)P , 


,7? 


77/ 


(85) 
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with 


CLL = - [L 2 + 4u p (u p - I 0 )] [L 2 + 4cn r (c J r + Sil 0 )] 0o(cos Opl, COS 6rl) 
CTL =pr(2u p - Z°) (—2sia; r - 1°) 4>i(cos9 pi , cos 64), 

L 2 

°TT = — (L 2 - 2 p 2 sin 2 - 2r 2 sin 2 9 ri ) 0 o (cos # pZ , cos 0 rJ ) 


2 

+4pV 2 0 2 (cOS0pj, COS &ri), 


where 


-2t r 


d(j) 


2n 


P 2 ( COS 0 pr ) (cos 0p r — COS 9 p i COS 9 r i) 


(f) n (cos 9 ph cos 9 rt ) = 

The explicit expressions we need are 

0o (x, y) = P 2 (x)P 2 (y), 

3 

0 iO , 2 /) = i ^?/ 0 2 - 1)(?/ 2 - 1 ), 

M x ,y) = TiX * 2 - i)(s / 2 - !)( 5 - 9a;2 - 9 y 2 + 2 izV). 

Id 


( 86 ) 

(87) 


( 88 ) 


Multiplying the resulting expression with Eq. ED we can read off the first 
contribution to 7i(p,r) from the box diagram: 

1 e 4 p r , 


= 


x 


dl 


3 -n F {u p ) [n B (w r - oj p ) + n F (w r )] 

57T' 3 j\fU p U r 

dl° (c TT |A^(L )| 2 + c LL |Af (L )| 2 + c TL Re A£(L)A£(L)) 


Wl 2 +AM 2 


x [«f(w p - /“) - n F (u r - I 0 )] 0(p - |L|)0(|/ + | -p)0(r - |L|) 6 (|( + | - r). 

(89) 

2. (si,s 2 ,s 3 ) = The cosines are cos 9 p i = Zpz, cos 9 r i = z~ t and the 

constraints are 


l 2 0 <l 2 , p > | 4 |, 

Therefore the second contribution reads 


r > |4|. 


(90) 


,(2) _ 1 e P r „/ 


KoxM = 


87t 3 Nf Up u x 


n' F (u p ) [n B (u r - Up) + n F (u r )} 


x 


dl 


dl° [c TT |A«(L )| 2 + c LL |Af (L )| 2 + c TL Re A£(L)A£ (L)) 


x [n F (n;p - /°) - n F (u r - Z 0 )] 0 (p - |Z+|) © 4 - |4l) • 


(91) 
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3. (si, s 2 ) S 3 ) = (+, —). In this case the cosines are cos 6 p i = z p i, cos 6 r i = z^ t , 
with the constraints 

ll < l\ P>|4l: r >\l-\. (92) 

The third contribution then reads 
1 4 

«Sc(p. r ) = ^jr—— n ' F M Mu r + w p ) + n F (w r )] 

o7T° iv f Up U r 

/ oo A 

dl I dl° [c TT |A^(L )| 2 + c LL |Af (L )| 2 + c TL Re A^(L)A^(L)j 
x [n F (l° + u r ) - n F (l° - Up)] 0(p- |Z+|)9(r — \l-\) . (93) 

Again, one can verify that Tibox is symmetric under exchange of p and r, which allows 
for the integral equation to be derived from the functional Q. 


5 Electrical conductivity 


The calculation of the electrical conductivity goes along the same lines as above. 
The coupling between two fermion fields and the external operator is simply 7 *. The 
correlator we need reads then 


Gjj(Q) = -q 2 f ^trS(P + Q)T i (P + Q,P)S(P)Y. 

Proceeding as we did for the shear viscosity, we introduce a scalar vertex as 
Reh A (p)r4p,cUp;p)-u A (p) = Reu A (-p)P(-a;p, -u p ; p)u A (-p) = 2 p 1 V(p), 
or through the equivalent expression 


V(p) = —Retr T*(p° + iO,p° — iO; p)(P + m)p l 
4 p 


p° = ±CJp 


We then find for the conductivity 


a = 


4 q 2 f N f 


P 2 Vf N f 


—n' F (up)x(p) = 


ip u p 


3n 2 


'0 


°° P 3 

dp — n' F (up)x{p)- 

UJp 


with 


x(p) 


P vip) 

r p 


(94) 

(95) 

(96) 

(97) 

(98) 
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The integral equation is actually simpler, since the contributions from the dia¬ 
grams with the box rungs cancel each other. This is a consequence of Furry’s theo¬ 
rem. It can be seen explicitly as follows. The integral equation for the conductivity 
has a similar form as Eq. (11)511 . but with P 2 (p • f) replaced by Pi(p ■ f) = p • r. Under 
the change of variables R —» —R mentioned below Eq. m, this factor is now odd 
and the second box diagram cancels exactly the first one. The integral equation reads 
therefore 

T\P + Q, P) = f + ^ ^7 V S{R + Q)T\R + Q, R)S(R)YD,„(R - P ). (99) 

Proceeding as before, we arrive at 

u P ?pX(p) =P + \ I [^(r 0 ) + n B {r° — o; p )] p • r^x(r )p(R)A(R, P), (100) 

z Jr r 

with the kernel 

M R , P ) = 2^-tr Kf + m W + m)7l PA R ~ P)- (101) 

Upon multiplication with 

2 

—4K), (102) 

Ldp 

the integral equation can be cast in the form of Eq. dZ2D, where P(p) is the same as 
above, but 

3 

S(p) = —n' F ( uj p ). (103) 

Ldp 

The kernel Tt receives only one contribution, 7Yii ne , which can be obtained from Eq. 
m by the replacement P 2 {z pr ) —> Pi(z pr ) = z pr . The electrical conductivity is given 
by 

Aq 2 f N f r / N 

o‘ =- bc = Xe*t\- (104) 

5.1 Ward identity 

The Ward identity relates the fermion-gauge boson vertex and the fermion self energy. 
We can explicitly verify that the integral equation which determines the fermion-gauge 
boson vertex is consistent with the Ward identity. We follow closely our analysis in 
the weakly coupled limit [1]. 
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The abelian Ward identity for the fermion-gauge boson vertex reads 

Q^{P + Q,P) — S~\P) - S~\P + Q ). (105) 

After performing the required analytic continuation and taking q = 0 this becomes 
g°rV + q° + iO+ : p° - i 0 + ; p) = q° 7 0 + S A (P) - 
I 11 terms of 

D(p) = kin g°^tr [T°(p° + q° + i0 + ,p° - A) + ; p)(f + M)\ 
the Ward identity in the pinching pole limit reads then 

®(p) = * r p- 

The integral equation for the zero’th component of the effective v< 

r°(P + Q,P) = 7 ” + A 1 I ‘ S ( R + Q) r °( R + e. R)S(Rh u D 

After performing the sums, the analytic continuation and taking appropriate limits 
as in Eq. (TTU7I) above, we find 

= J R ^ r ° _ _ nF ^ 

x {2 p T (R — P ) [r°p° — m 2 — k • pk • r] + Pl(R — P) \p°r° + m 2 + p ■ r] j , (110) 

where p° = ±u; p . It is now straightforward to see that the solution of the previous 
integral equation is precisely given by the Ward identity, Eq. (HUD, since after substi¬ 
tuting D(r) = fT r inside the integral we obtain exactly the expression for the thermal 
width Eq. (ED. We End therefore that the diagrammatic formulation is consistent 
with the Ward identity. 

We note here that the analysis of the Ward identity is easier in the large Nf 
limit than in the weakly coupled limit in the leading logarithmic approximation |5]. 
The reason is that the soft-fermion contribution in the leading-log approximation, 
which leads to an additional diagram in the integral equation (see Ref. (5) for further 
details), is subleading in the large Nf limit and does not have to be considered. 


P + Q )• (106) 


p° = ±LJ p 


(107) 


ertex is 
pu(R-P). 


(108) 


(109) 
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6 Color factors 


The results for large Nj QCD can be obtained from the previous analysis, provided 
that the appropriate color factors and electric charges of the quarks are inserted. The 
relevant group factors are 

N 2 — 1 1 

Cp = 2N ’ Tf = 2 ’ dp = N °' ( m ) 

In terms of the effective coupling [p 

9ls = 9 2 T F , (112) 


where g is the QCD coupling rescaled with 1 / yJTTf, so that it remains finite in the 
large Nf limit (cf. Eq. (|2|) ). the QCD expressions follow from the abelian results with 
the replacements 


quark self energy/thermal width: 
gluon self energy: 

“line” piece in the integral equation: 
“box” piece in the integral equation: 


2 Cf 

d eff rj~i 

1 F 

(113) 

9eS 

(114) 

2 Cp 

9eS rp 

1 F 

(115) 

4 Cp 

9e ff rp 
±F 

(116) 


If we now define x as (compare with Eqs. (1641981) ) 

P 2 V(p) C F 


xip) = 
xip) = 


(Up T p Tp 
p V(p) Cp 
Up T p Tp 


(shear viscosity), 
(electrical conductivity), 


(117) 

(118) 


the integral equations are exactly the same as in QED, with no explicit color factors. 
The transport coefficients are then given by 


and 


V = 


a = 


dpTp 2 Nf 


dpTp 


Cp 3n 2 

where qj is the electric charge of the quarks. 


f°° p 4 

/ dp — n' F (u p )x(p), 

Jo ^P 

(119) 

- [ dp n' F (u p )x(p), 

Jo ^ p 

(120) 
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7 Variational solution 


In order to obtain the shear viscosity and electrical conductivity for general values of 
the effective coupling constant and mass parameter, we solve the problem of extrern- 
izing the functionals Q variationally. Following Refs. m 131 X 3 , we expand x(p) hi a 
finite set of suitably chosen basis functions (j)^ m \p): 

N v ar 

X(p) = N fYl ( 121 ) 

m=l 

where we factored out an explicit factor of Nf, so that the integrals below are Nf 
independent. With this Ansatz, the functional Q reads 


with 


«M = iV/E' 


*5 m „ ^ ( •Ft r 


mn i ' '-mn 


H n 


( 122 ) 


c 


dp S(p)d, {m \p), 


r m „ = n, dp T(p)^ m] (p)d> {n) (p), 

Jo 

/•OO /*oo 

’Hmn = N f dp / dr’H(p,r)</> (m> (p)</> l - n ' l (r). 


(123) 


Extremizing the functional with respect to the variational parameters a m gives the 
solution 

Gm = E(^-^)-A ( 124 ) 


so that the shear viscosity and electrical conductivity are given by 

dpTp 


dpTp 2 Nf 


a = 


Cp 37r 2 


E 5 " 


(125) 


with the corresponding values of S m and a m for each transport coefficient. S m is a 
1-dimensional integral, T mn and are 3-dimensional integrals and for the shear 
viscosity is a 4-dimensional integral. These integrals are done with numerical 
quadrature. The uncertainty due to the numerical integration is estimated to be on 
the percent level or less. 
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We now discuss the choice of trial functions. We work with the set 


(126) 

(127) 


0 (m) (p) = (p/T) 2 ft m \x) for rj, 

0( m )(p) — p/T 2 ft m \x) for a. 

where x = p/\/T 2 + mT. The prefactors of p are motivated by an approximate 
analytical solution of the integral equation for massless fermions in the leading log 
approximation for asymptotically large momentum (see Ref. m for a similar analysis 
for Nf = 1). The dimensionless argument x is chosen such that it represents the 
magnitude of the typical momentum, which is p ~ T for small fermion masses and 
p ~ a JmT for larger masses. For the dimensionless functions (p' rn ' ! (x) we use 



1 

(1 + x) m ~ 1 


m— 1 

£(-*)* 


k =0 


(128) 


for not too large fermion masses, m/T < 5. This set was also used in Ref. [EEJ for 
the calculation of the shear viscosity in the O(N) model. For larger mass, however, 
we find that this choice of trial functions leads to slow convergence in the number of 
trial functions. In this case we use Laguerre polynomials, <p m \x ) = L m _i(x). The 
results shown in the figures are obtained with a set of dimension 4, and the effect of 
using this truncated basis set is smaller than the width of the lines. 

The dependence of the shear viscosity on the fermion mass and the effective gauge 
coupling is shown in Figs. [Till and EH while the results for the electrical conductivity 
are shown in Figs. H21andlT31 The results are normalized with 


Vo 


d F T F Nj 
Cf g 4 eS 


a 0 = 


d F T F 
C F g 4 s 


(129) 


For QED, g eS = e and d F = T F = C F = 1, while for QCD, g^ s = g 2 T F . We remind 
that both e and g remain finite in the large Nf limit and that for the purpose of 
presentation we have chosen p = hdr = ire~ lE T. 

We notice that the general behavior of these transport coefficients is to decrease 
with increasing mass except for small values of the coupling constant, where a slight 
increase for small masses is observed. After rescaling with g o resp. oo, the remaining 
dependence on the effective coupling constant is quite strong, much stronger than 
the subleading dependence on the quartic coupling in the case of the shear viscosity 
in the O(N) model [I2J. In the limit of vanishing fermionic mass, our results agree 
with those obtained in Ref. 2] using kinetic theory. In the opposite limit of very 









Figure 10: Shear viscosity vs. the fermion mass m for various values of g e g. 



Figure 11: Shear viscosity vs. the effective coupling constant g e g for various values of 
the fermion mass m. 
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Figure 12: Electrical conductivity vs. the fermion mass m for various values of g e g. 



Figure 13: Electrical conductivity vs. the effective coupling constant g e g for various 
values of the fermion mass m. 
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large mass, we provide in Appendix El some parametric estimates of these transport 
coefficients in the leading logarithmic approximation, see Eq. (HMD, which corroborate 
the behavior shown in the plots. The difference between the mass dependence of the 
shear viscosity in the O(N) model [T3] and in the gauge theory is due to the different 
mass and momentum dependence of the (transport) cross section. 

In our complete leading order calculation in the large Nf limit, the presence of 
the Landau pole means that, for given fermion mass, there is an upper bound on 
the possible values of the coupling constant, arising from the requirement that all 
physical scales lie well below the Landau scale Al- For large masses m/T 1, the 
typical momentum scale is not set by the temperature, but instead by p rsj \JvnT S> T. 
Thus, the upper bound on the coupling constant decreases when the mass is increased. 
This implies that going to asymptotically large mass requires a restriction to the weak 
coupling limit. 

8 Conclusions 

We have presented a diagrammatic calculation of the shear viscosity and the electrical 
conductivity at leading order in the large Nf expansion of QED and QCD for massive 
fermions. 

The 2PI effective action at next-to-leading order provides in a straightforward 
manner appropriate integral equations which sum all the required diagrams to obtain 
these transport coefficients at leading order. We proved that these equations are 
gauge fixing independent and consistent with the Ward identity. This explicitly shows 
that in a fully self-consistent calculation of these transport coefficients at leading 
order in the 2PI framework, potential non gauge invariant contributions would be 
suppressed by powers of Nf. This suggests that in self-consistent applications of the 
2PI effective action in gauge theories, e.g. as in far-from-equilibrium applications, 
potential problems related to gauge invariance and Ward identities would be small 
for sufficiently large Nf. 

Our results show a nontrivial dependence of the shear viscosity and electrical 
conductivity on the mass of the fermions and the effective gauge coupling. We found 
that for small values of the coupling constant they increase slightly with increasing 
mass. For larger values of the fermion mass, both the shear viscosity and electrical 
conductivity decrease. It would therefore be interesting to extend current calculations 
of transport coefficients to include different massive fermion flavors. We also found 
that after taking out the expected 1/a 2 dependence, a strong dependence on the 
gauge coupling remains. When nonperturbative results obtained from lattice QCD 
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simulations eh are compared with perturbative ultrarelativistic expressions, these 
findings should be kept in mind. 
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A 2PI effective action 


In this appendix we summarize some useful exact relations derived from the 2PI 
effective action. We consider only bilocal sources, such that the path integral is 

Z{K,H\=e iW[K ’ H] = j V$VipV(t) e <^^ ] + 5 ^y‘P+^rb i ') _ ( 130 ) 

We denote Bose fields collectively with 0* and fermion fields with 0 a . Indices indicate 
space-time as well as internal indices, and integration and summation over repeated 
indices is understood, e.g., 


H ba '^ a '^ b = I H^ a (y, x)^ a (x)^(y) = I tr H(y, x)ip(x)ip(y). 


’xy 


' xy 


The 2PI effective action follows from the Legendre transform 


where 


r [G, S} = W[K , H] - -G^K l3 - H ba S ab , 


SW 1 .. 1 

GU= _ (W >, 


L u 

5W 


= s ab = (T c r^ b ), 




SH ba 

The effective action is written as 


JL = _!*.. 

8Gv 2 lv 


5T 
6S Va 


= -H nh . 


r[G, S] = -Tr In G _1 + -Tr G^G - G 0 ) - iTr In 5” 1 


(131) 


(132) 


(133) 

(134) 
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zTySo 1 (S-S 0 ) + T 2 [G, S], 

(135) 






where 


<G -i = JjL 

013 SfiSft 


(j)=ljj='ljj=0 

Varying this effective action once yields 

^ = -l (G- 


j Q —1 — 
ab ~ 


8 2 S 


S'i/j a 8 , ijj b 


(j)=ijj=tl )=0 


SG* 


1-1 

r ij 


G olj + n *j), 


sr 


with 


„ „. 8T 2 

Uij 1 JG h 


8S ba 

^ab ^ 


~ * (^aft ‘S'oab + J 

£r 2 




Varying one more time results in the relations 


- 7 { G ik G jl + G u 1(G jk ~ ^ij-M ) , 


8G ij 8G kl 4 

5 2 r ! 

figbafigdc — * (~^ad &cb ~ ^-ab;cd) j 

5 2 r i 


8G ij 8S ba 


= -Ain- 


ij\ab ? 


where we dehned 

. 5 2 r 2 
— 4i- 


' 8G ij 5G kl ’ 

From Eqs. cm SUP we find 

5 2 IF 


Aij-ab 2z 


< 5 2 r 2 


8G^8S ba ' 


A a b\cd f 


< 5 2 r 2 


5S ba 8S dc ' 


S K,,mri iGfM+G ' kG,l+G ' ,G,k) ’ 


8 2 W 


8H ba 5H dc 

5 2 W 


= i (Gf’ cd - S cb S ad ), 


SKijSHba ~ 2 c ’ 

where G l J' ,H etc. are the usual connected 4-point functions. It is convenient 
vertex functions dehned by truncating legs 

GfAi’ = G*G ij, G kk 'G 11 ' rg, ;fc , z ,., 

Qij-,ab _ G™' S aa ' m 

f'iab\cd\ naa' nb'b qcc' nd'dr^i 4) 

Lt c — O OOO L a'b'-, C 'd'- 


(136) 

(137) 

(138) 

(139) 

(140) 

(141) 

to use 

(142) 
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From the following identities 


n JJtij 

5HJ 


rdcc — ab 

°a°h ~ 


SH, ’ 


°=fy (i«) 


‘6a V fide u x v zj 

we arrive at four coupled integral equations for the 4-point vertex functions: 


p(4) _ A I * A A QCC' Qfl'flrW 

L ij;ab ~ IV ij\ab ~r Lt 1 k'l';ab lv ij',dc>J 1 c'd';a6i 

rJtt = W + - A 0 j 1/e s“'s / ' / r^}, ;c j. 

p( 4 ) _ * , }_ A p'ffcfc'WZ'pP) _ » qcc' cd'rfpl 4 ) 

1 ab;ij ~ 1V ab;tj T ^ab-,kl^ Lt 1 k'l';ij Iy ab;dcJ *-> 1 c 'd';rp 


ittV«'r(4) 


cc' Qcfc(p(4) 


p(4) a I " A rimin'pmn'pW a qcc' Qd'dpi 4 ) /'i/l/l'l 

1 ij;kl ■'Vjj.y + ^Ajj. mn (jr Or *- m ' n '-kl ^Hj^dc^ *-> *- c 'd'\ij- l 144 ! 

In the main text we employ these equations using the 1/Nf expansion of the 2PI 
effective action at NLO. 




44) 


cc' Qd'dp(4) 


B Parametric estimates 


In this appendix we discuss parametric estimates in the leading logarithmic approx¬ 
imation, in the zero and large fermion mass limit, using a standard kinetic theory 
discussion (see e.g. |22j). 

It follows from the hydrodynamical definitions 0 that the shear viscosity and 
the electrical conductivity are related to a diffusion coefficient as 


D 


T/ 

{s + vy 



(145) 


where £ is the energy density, V the pressure, and S the charge susceptibility. For 
parametric estimates, this diffusion constant can be taken to be the same, since similar 
processes determine the transport of energy momentum and charge in the large Nf 
limit, namely large angle scattering between fermions. The diffusion constant can be 
estimated using a random walk model [2H ES, as D = £ m fV, where v is the average 
speed and ~ l/h<7 tr the mean free path, with n the mean density and a tI the 
transport cross section Ea 


&tT — 


d£l{ 1 


cos 6) 


dcr 

dQ 


(146) 
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If we keep only the most divergent term in the differential cross section, 
dir a 2 1 1 da a 2 m 2 1 

dn~Wfs¥ (p>>m) ' dn~N]V¥ (p<<m) ' 

where a = e 2 /47r, and we recall that the gauge coupling e 2 is rescaled with 1 /Nf, the 
transport cross section reads 



O'tr ~ 


N 2 T 2 



a 


3-1 


In 0~ 
min ‘ 


(148) 


Here we used that in the relativistic case the typical momentum p ~ y/s ~ T while 
in the nonrelativistic case p ~ mv ~ y/Tm. Eq. (|148j) holds for both p m and 
p -C m. The divergence at small angles is, to leading log accuracy, cut off by Debye 
screening m 

For light fermions (T m) we use that 
{£) ~ (V) ~ NfT\ n ~ NfT 3 , S ~ N f q 2 f T 2 , v ~ 1, (149) 


as well as that # min ~ q/p ~ mp/T ~ e, where q is the exchanged momentum. This 
leads to the well-known parametric estimates mm 


V 


Ni 


T 3 


cr 


In 1/a ’ 


cr 




T 


a 2 In 1 /a 


(150) 


In the case of heavy fermions ( m/T 1) in the regime where scattering can be 
treated classically ( a 2 m/T 1) ;22], the exchanged momentum q can be estimated 
as the product of the force a/r 2 D and the transit time r D /v for a passage at impact 
parameter ro, the typical Debye distance [22 . This yields q ~ a/(r£>v). We find 
therefore that # min ~ q/p ~ a/(roT). The inverse Debye mass td = 1 /mo is deter¬ 
mined from m 2 D = — n°°(p° = 0, p —> 0). For large m>Twe hnd an exponentially 
small Debye mass, 


m 2 D ~ am 2 



e -m/T [J + Q ( T/m )] , 


(151) 


such that 


log# 


-i 


m 3, T 

TT + 7 log W - 

21 4 arm 


+ 0(T/m). 


(152) 


S ~ q)n/T, 


Combining this with 

(£) ~ mn, 


{V) ~ Tn, 
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V r ^ 1 Vfj m, 


(153) 









yields 


V 


T 3 

N f— 

J cr 


T \2 

m J 


<7 ~ 




(154) 


As indicated above, these expressions are valid in the leading log approximation 
log T/amr, 1, the large mass limit m/T 1, as well as the classical scattering 
limit a 2 m/T 1. 
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